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Abstract 

We prove the existence and uniqueness of a strong solution of a stochastic differential 
equation with normal reflection representing the random motion of finitely many globules. 
Each globule is a sphere with time-dependent random radius and a center moving according 
to a diffusion process. The spheres are hard, hence non-intersecting, which induces in the 
equation a reflection term with a local (collision-)time. A smooth interaction is considered 
too and, in the particular case of a gradient system, the reversible measure of the dynamics 
is given. In the proofs, we analyze geometrical properties of the boundary of the set in 
which the process takes its values, in particular the so-called Uniform Exterior Sphere and 
Uniform Normal Cone properties. These techniques extend to other hard core models of 
objects with a time-dependent random characteristic: we present here an application to 
the random motion of a chain- like molecule. 

AMS 2000 subject classification: 60K35, 60J55, 60H10. 

Keywords: Stochastic Differential Equation, hard core interaction, reversible measure, 
normal reflection, local time, Brownian globule. 



1 Introduction 

Since the pioneering work of Skorokhod [TU], many authors have investigated the question 
of the existence and uniqueness of a solution for reflected stochastic differential equations 
in a domain. It has first been solved for half-spaces, then for convex domains (see 
Lions and Sznitman [5] proved the existence of a solution in so-called admissible sets and 
Saisho [8] extended these results to domains satisfying only the Uniform Exterior Sphere 
and the Uniform Normal Cone conditions (see definitions in Section \'d.l\) . These results 
were applied to prove the existence and uniqueness of Brownian dynamics for hard spheres 
in [6], [2], [3], or for systems of mutually reflecting molecules (see [9]). 

We are interested here in dynamics of finitely many objects having not only a random 
position but also another random time-dependent geometrical characteristic like the radius 
for spheres or the length of the bonds in a molecule. We will prove the existence and 
uniqueness of random dynamics for two elaborated models, using methods which are re- 
finements of Saisho's techniques, analyzing fine geometrical properties of the boundary of 
the domain on which the motion is reflected. 

More precisely : 
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We first introduce a globules model representing a finite system of non-intersecting 
spheres which centers undergo diffusions and which radii vary according to other diffusions 
constrained to stay between a maximum and a minimum value. The spheres might be cells, 
or particles, or soap bubbles floating on the surface of a liquid (2-dimensional motion) or 
in the air (3-dimensional motion). The behavior of the globules is quite intuitive : two 
globules collide when the distance between their random centers is equal to the sum of 
their random radii, and the collision has as effect that their centers move away from one 
another and their sizes decrease. The associated stochastic differential equation (£ g ) (see 
Section l2TTj) includes several reflection terms, each of them corresponding to a constraint in 
the definition of the set of allowed globules configurations : constraints of non-intersection 
between each couple of spheres, constraints on the radii to stay between fixed minimum 
and maximum values. We proves that this equation has a unique strong solution and give 
in some special time- reversible initial distribution (theorems 12.21 and 12. 4p . 

We also consider a model for linear molecules, such as alkanes (carbon chains) or poly- 
mers : each atom moves like a diffusion, the lengths of the bonds between neighbour atoms 
vary between a minimum and maximum value which evolve according to a reflected diffu- 
sion. This corresponds to a SDE (E c ) reflected on the boundary of the set of all allowed 
chains. Here also, we prove the existence and uniqueness of the solution of (£ c ) (see theorem 
12. 5p with similar methods as in the globule case. 

The rest of the paper is organized as follows : In Section [3j we present a new general 
criterion for a multiple constraint domain to satisfy the Uniform Exterior Sphere and 
Uniform Normal Cone conditions (see proposition I3.4p . These geometrical assumptions 
on the boundary induce existence and uniqueness of the reflected diffusion on this domain. 
We also obtain a disintegration of the local time along the different constraint directions . 
Section H] is devoted to the proofs of the theorems announced in Section [2j 

For the sake of shortness, we restricted ourselves to these two examples, though dy- 
namics for Brownian systems evolving under multiple constraints may be found in other 
situations (the results in section [3] are given in a general frame for easier adaptation to 
other examples). 

2 Two hard core models 
2.1 Globules model 

We want to construct a model for interacting globules. Each globule is spherical with 
random radius oscillating between a minimum and a maximum value. Its center is a point 
in M. d , d > 2. The number n of globules is fixed. Globules configurations will be denoted 
by 

x = (x\, xi, . . . , x n , x n ) with xi,...,x n GM d and xi,...,x n EM. 

where is the center of the i th globule and Xj is its radius. An allowed globules configuration 
is a configuration x satisfying 

Vz r_ < Xi < r + and Mi ^ j \xi — Xj\ > Xi + Xj 

So, in an allowed configuration, spheres do not intersect and their radii are bounded from 
below by the minimum value r_ > and bounded from above by the maximum value 
r + > r_. In this paper, the symbol | • | denotes the Euclidean norm on W 1 or R( d+1 ) n (or 
some other Euclidean space, depending on the context). 
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Let A g be the set of allowed globules configurations : 

A g = {x G IR (d+1)n , Vz r_ < < r + and Vz 7^ j |x, - Xj| > X; + x,} 

The random motion of reflecting spheres with fluctuating radii is represented by the fol- 
lowing stochastic differential equation : 

X t (t)=X t (0)+ f <Ti{X(8))dW t {8) + f\{X(a))da + Y i f * { *\~ f^M s) 

JO Jo X ii S ) + X ji S ) 

rt rt n 

X l (t)=X l (0) + / #i(X(s))dWi(s) + / %(lt(s))ds-Y^L i3 (t)-Lj(t) + LT(t) 
{ Jo Jo j=1 



In this equation, X(s) is the vector (Xj(s), Xj(s))i<j<„. The Wj's are independents Re- 
valued Brownian motions and the Wj's are independents one-dimensional Brownian mo- 
tions, also independent from the W^s. The diffusion coefficients <7j and cfj, and the drift 
coefficients 6, and 5, are functions defined on A g , with values in the d x d matrices for 
(7j, values in K d for 63, and values in R for <jj and 6j. To make things simpler with the 
summation indices, we let L# = 0. 

A solution of equation (£ 9 ) is a continuous *A 9 -valued process {X(£),i > 0} satisfying 
equation (S g ) for some family of local times Ly, Lf, L~ such that for each i,j : 



Lij — Lji, Lij(t) — / l|x i ( s )-x j (s)|=x i (s)+x J (s) dLij(s) , 

t t 
Lt(t) = [ l Us)=T+ dLf{s) and Lr(t) = ! 1^ dL.r( s ) 
JO •/ 



Remark 2.1 : Here, and in the sequel, the expression local time stands for : non- 
decreasing adapted continuous process, which starts from and has bounded variations on 
each finite interval. 

Equation (E g ) has an intuitive meaning : 

• the positions and radii of the spheres are Brownian; 

• when a globule becomes too big, it's deflated : Xi decreases by —dLf when Xi = r + ; 

• when a globule becomes too small, it's inflated : Xi increases by +dL~ when Xi = r_; 

• when two globules bump into each other (i.e. \Xi — Xj\ = Xi + Xj), they are deflated 
and move away from each other : Xi decreases by —dL^ and Xi is given an impulsion 
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Xi—Xj 



in the direction ^\^ 3 with an amplitude dLij 



In the case of an hard core interaction between spheres with a fixed radius, the existence 
of solutions for the corresponding SDE has been proved in [6j. However, the condition 
\ x i ~ x j\ > %i + Xj is not equivalent to | (xi, Xj) — (xj ,Xj)\ > c for some real number c. Hence 
the above model is not a classical hard sphere model in IR d or IR d+1 . 

Theorem 2.2 Assume that the diffusion coefficients o~i and <Xj and the drift coefficients bi 
and hi are bounded and Lipschitz continuous on A g (for 1 < i < n). Then equation (E g ) 
has a unique strong solution. 
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Remark 2.3 : "Strong uniqueness of the solution" here stands for strong uniqueness 
(in the sense of [I] chap. IV def.1.6) of the process X, and, as a consequence, of the reflection 
term. This does not imply strong uniqueness for the local times L^, Lf, L~ unless several 
collisions at the same time with linearly dependant collision directions does not occurs a.s. 
(see the proof of corollary 13.61 for details). 

The first part of the next theorem is a corollary of the previous one. The second part 
describes the equilibrium states of systems of interacting globules. Here, and through this 
paper, <ix denotes the Lebesgue measure. 

Theorem 2.4 The diffusion representing the motion of n globules submitted to a smooth 
interaction (p exists as soon as the interaction potential if is an even C 2 function on M. d 
with bounded derivatives. It is the unique strong solution of the equation : 



I ft n n r 

X l (t)=X l (0) + W t (t)-- / Y,^^Ms)-x,(s))ds + J2 \ 

ZJ ° 3=1 1=1 J ° 



X z (s)+X,(s) 



dLi^s) 



x { (t) = 1,(0) + w t {t) - 5>«(f) - + L T(t) 

3=1 

with X Ag-valued continuous process and (L^, Lf, Lj)o<i,j<n satisfying conditions (£') 

Moreover, if Z = e~^-' 1 ^ i< i^ n ^ Xi ~ x ^d'x < +oo then the Probability measure /x defined 
by d/x(x) = ^l_4 s (x)e~ Si<i<j< n ^( X4_x j)(ix is time-reversible for this diffusion, i.e. (X(T — 
t))te[o,T] has the same distribution as (K(t)) t e[o,T] for each positive T when X(0) ~ ji. 



These theorems are proved in section 14.11 The proofs rely on previous results from 
Y.Saisho and H. Tanaka (see section I3TTI) and on an inheritance criterion for geometrical 
properties which is given in section 13.21 



2.2 Linear molecule model 

Another example of hard core interaction between particles with another spatial charac- 
teristic beside position is the following simple model for a linear molecule. In this model, 
we study chains of particles having a fixed number of links with variable length. More 
precisely, a configuration is a vector 

x = (xi, . . . , x n , x x + ) with ii,...,i n 6 M d and x_,x + GM 

Xi and Xi+x are the ends of the i th link and x_ > (resp. x + > x_) is the minimum (resp. 
maximum) allowed length of the links for the chain. The number n of particles in the chain 
is at least equal to 2, they are moving in M d , d > 2. So the set of allowed configurations 
is : 

A c = {x e R dn+2 , r_ < x_ < x + < r + and Vz G {1, . . . , n - 1} x_ < \xi - x i+1 \ < x + ] 
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We want to construct a model for the random motion of such chains, as the ^4 c -valued 
solution of the following stochastic differential equation : 



X t (t)=X l (0) + [ ai(X(s))dWi(s) + f b t (X(s))ds+ [ Xi Xl+l (s)dL; (s) 
Jo Jo Jo X- 

+ f^J^{s)dL-_ x {s)- f^J^{s)dLt{s)- f ^J^{s)dLUs) 
Jo yi- Jo X + Jo X + 

X_(t)=X_(0)+ / <r_(X(s))dW_(s) + / L(X(s))ds-Vl:(i) + L_(<)-L = (t) 
</o Jo i=1 

X + (t)=X40)+ / a + (X(s))dW + (s) + / 6 + (X( S ))d S + ^L+(t)-L + (t) + L = (t) 
./o Jo i=1 

As in the previous model, the W^s are independents l^-valued Brownian motions and the 
W- and W + are independents one-dimensional Brownian motions, also independent from 
the Wj's. A solution of equation (£ c ) is an *A c -valued continuous process {X(t),t > 0} 
satisfying the equation for some family of local times , Lf, L_, L + , L = such that for 
each t G M + and each 2 : 

L i (*) = ^ I |X l (s)-X l+1 ( S )|=X_( S ) C?L i ( S ) > = ^ I |X l ( S )-X l+1 ( S )|=X + ( S )^ + (- S ) 

L-(t)= [ l^ {s)=r _dL4s), L+(t)= f l x+{s)=r+ dL + (s) and L = (t) = f l x _ {s)=x+{s) dL = (s) 
Jo Jo Jo 

Equation (S c ) looks more complicated than (S g ) because it contains a larger number of local 
times, but it is very simple on an intuitive level : both ends of each link are Brownian, 
links that are too short (\Xi — Xi+\\ = X_) tend to become longer (reflection term in 
direction Xi —Xi + \ for Xi and in the opposite direction for Xj + i) and also tend to diminish 
the lower bound X_ (negative reflection term in the equation of X_). Symmetrically, 
links that are too large (|Xj — X i+ i\ = X + ) are both becoming shorter (reflection term in 
direction X i+1 — Xi for Xi) and enlarging the lower bound X + (positive reflection term 
in its equation). Moreover, X_ increases (by L_) if it reaches its lower limit r_ and X + 
decreases (by L + ) if it reaches its lower limit r+. And the lower bound decreases and the 
upper bound increases (by L=) when they are equal, so as to fulfill the condition X- < x + . 

Theorem 2.5 If the a^s, c\_ and a + and the bi's, b_ and b + are bounded and Lipschitz 
continuous on A c , then equation (S c ) has a unique strong solution. 

Moreover, assume that the Oi 's, cr_ and a + are equal to the identity matrix, 6_ and b + 
vanish, and 6j(x) = — | X^j=i ^f( x i ~ x j) f or some even C 2 function ip on M> d with bounded 
derivatives, satisfying Z = e~^- ,1 ^ i <^ n ' p ^ Xi '' Xj 'dyc < +00. Then the solution with initial 
distribution dfi(x.) = ^l_4 c (x)e~ ^^ !< ^ n ^^"^dx is time-reversible. 



See section 14.21 for the proof of this theorem. 
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3 Geometrical criteria for the existence of reflected 
dynamics 



3.1 Uniform Exterior Sphere and Uniform Normal Cone proper- 
ties 

In order to solve the previous stochastic differential equations, we will use theorems of 
Saisho and Tanaka extending some previous results of Lions and Sznitman [5] . 

To begin with, we need geometrical conditions on subset boundaries. The subsets we 
are interested in are sets of allowed configurations. But for the time being, we just consider 
any subset T> in R m (m > 2) which is the closure of an open connected set with non-zero 
(possibly infinite) volume, dx is the Lebesgue measure on R m , | • | denotes the Euclidean 
norm as before, and we set 

5 m = {xer, |x| = l} 
We define the set of all (inward) normal vectors at point x on the boundary &D as : 

Af? ={J where N? >a = {n G S m , 5(x - an, a) n V = 0} 

Here B(pc,r) is the open ball with radius r and center x. Note that 

S(x - an, a) n V = Vy G V (y - x).n + — |y - x| 2 > 

2a 



Definition 3.1 If there exists a constant a > such that N"^ = ^ for each 

x G &D, we say that T> has the Uniform Exterior Sphere property (with constant a) 
and we write : V G UES(a). 

UES(a) means that a sphere of radius a rolling on the outside of T> can touch each point of 
&D. This property is weaker than the convexity property (which corresponds to UES '(oo)) 
but still ensures the existence of a local projection function similar to the projection on 
convex sets. 

Definition 3.2 We say that T> has the Uniform Normal Cone property with constants 
(3, 5, and we write V G UNC(/3,8), if for some (3 G [0, 1[ and 5 > 0, for each x G &D, 
there exists l x G S m such that for every y G dT> 

| y - x | < 5 =^ Vn G n.l x > Vl-/? 2 

Let us consider the reflected stochastic differential equation : 

X(t)=X(0)+ [ (r(X(i))dW(s) + / b(X(s))ds+ f n(s)dL(s) (1) 
Jo Jo Jo 

where (W(t)) t is a m-dimensional Brownian motion. A solution of ([T]) is a (X, n, L) with 
X a P-valued adapted continuous process, n(s) G A/^/ s \ when X(s) G dV and L a local 
time such that 

L(-) = / lx(s)edvd~L(s) 
Jo 

The following result is a consequence of [B] and [7] : 
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Theorem 3.3 Assume T> G UES(a) and D e UNC(/3, S) for some positive a, S and some 
j3 6 [0, 1[. If a : X> — > R m and b : T> — > W 71 are bounded Lipschitz continuous functions, 

then (CP has a unique strong solution. Moreover, if er is the identity matrix and h = — -V$ 

with $ a C 2 function on IR m with bounded derivatives satisfying Z = f v e~^^dx < +oo, 
then the solution with initial distribution d/i(x) = -|lx)(x)e _ * < - x - ) (ix is time-reversible. 

Proof of theorem 13.31 

The existence of a unique strong solution of ([T]) is proved in [5] theorem 5.1. 

Let us consider the special case of the gradient system : b = — |V$ and cr identity matrix. 
For fixed T > and x 6 X>, let P x be the distribution of the solution (X, n, L) of 
starting from x. The process (W(t)) t ^[o,T] defined as 

W(t) = X(t) -x+~ / V$(X(a))ds- / n(s)dL(s) 
2 Jo Jo 

is a Brownian motion with respect to the probability measure P x . Since V$ is bounded, 

1 /•* 

thanks to Girsanov theorem, the process W(t) = W(t) / V$(X(s))ds is a Brownian 

2 Jo 

motion with respect to the probability measure P x defined by — — = where 

aPx 

M t = expQjT V$(X(s)).dW(s) - ijT |V$(X( S ))| 2 ^ 
From theorem 1 of [7], it is known that Lebesgue measure on X> is time-reversible for the 



solution of X(i) = X(0) + W(t) + / n(s)c?L(s), that is, the measure Pd x = / P x <ix is 

Jo Jv 
invariant under time reversal on [0; T] for any positive T. Using Ito's formula to compute : 



f V$(X(«)).(iW(g) = *(X(T))-*(X(0))+J / |V$(X(s))| 2 -A$(X(s))(is- / V$(X(s)).n(s)dL(s 
Jo 2 j Jo 



we notice that the density of the probability measure P M = / P x //(dx) with respect to 
the measure Pi x is equal to : 
-*(X(o)) /$ ( x(0))-*(X(T)) , ri A ^^ l™ v ,^ l2 ^ , 1 r. 



■exp( $(X(Q)) 2 HMT)) +1 ^(X( S ))-i|Vd>(X(,))| 2 d S + ^ V*(X(*)).nOOdL( a )) 



This expression does not change under time reversal, i.e. if (X, n, L) is replaced by 
(X(T — ■), n(T — ■), L(T) — L(T — ■). Thus the time reversibility of P^ x implies the time 
reversibility of P M . ■ 



3.2 The special case of multiple constraints 

The sets we are interested in are sets of configurations satisfying multiple constraints. They 
are intersections of several sets, each of them defined by a single constraint. So we need 
a sufficient condition on the XV s for Uniform Exterior Sphere and Uniform Normal Cone 
properties to hold on X> = flf =1 XV 

In [H], Saisho gives a sufficient condition for the intersection X> = nf =1 X , i to inherit both 
properties when each set X>j has these properties. However, in order to prove the existence 
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of a solution in the case of our globules model, we have to consider the intersection of the 
sets : 

Vij = {x G R dn+n , \xi -Xj\>Xi + xj} 

Uniform Exterior Sphere property does not hold for X^-'s, so Saisho's inheritance criterion 
does not work here. However, T>ij satisfies an Exterior Sphere condition restricted to A g in 
some sense, and we shall check that this is enough. Similarly, the Uniform Normal Cone 
property does not hold for V^, but a restricted version holds, and proves sufficient for our 
needs. 

So we present here a UES and UNC criterion with weaker assumptions. We keep 
Saisho's smoothness assumption, which holds for many interesting models, and restrict the 
UES and UNC assumptions on the TVs to the set T>. Instead of the conditions (Bq) and 
(C) in [9], we introduce the compatibility assumption (iv) which is more convenient because 
it is easier to check a property at each point x G dT> for finitely many vectors normal to 
the &DiS, than on a neighborhood of each x for the whole (infinite) set of vectors normal 
to &D. 

The proof of this criterion is postponed to the end of this section. 

Proposition 3.4 (Inheritance criterion for UES and UNC conditions) 

v 

For V C M. m equal to the intersection V = f^T^i, assume that : 

i=l 

(i) The sets V>i are closures of domains with non-zero volumes and boundaries at least C 2 
in T> : this implies the existence of a unique unit normal vector rij(x) at each point 
x eVDdVi. 

(ii) Each set T>i has the Uniform Exterior Sphere property restricted to T>, i.e. 
3«i > 0, VxG£>n<9£>; B(x- otiDiipc),^) HV { = 

(Hi) Each set T>i has the Uniform Normal Cone property restricted to T>, i.e. 

for some fa G [0, 1[ and 5i > and for each x G T> fl dT>i there is a unit vector 1^. s.t. 
VyeVn&DiHBfaSi) n 4 (y)4 > y^T^f 



(iv) (compatibility assumption) There exists (3$ > ^J2 maxi<j< p /3j satisfying 
Vx G &D, 31° G S m , \/i s.t. x G dVi l°.n,(x) > fa ~ 

Under the above assumptions, T> G UES(a) andT> G UNC (fa 5) hold with a = fa mini<j< p ai, 
5 = min 1 <j< p 5j/2 and (3 = \J\ — (fa — 2max 1 <j< p /5j//5 ) 2 . Moreover, the vectors normal 
to the boundary dV are convex combinations of the vectors normal to the boundaries dV>i : 

Vx G dV JV£ = I n G S m 1 n = ^ qn^x) with each a > I 

Remark 3.5 : Thanks to the compatibility assumption (iv), n = Xlac 9x c * n *( x ) with 
non-negatives q's and |n| = 1 implies that 

E « < E ^ )A ° X = — < 1 

so that the last equality in proposition 13.41 can be rewritten as 

Af? = < n G S m , n = E c * n *( x ) with V * Q > and E Q - }f f 



8 



Corollary 3.6 of theorem 13.31 

v 

If T> = O T>i satisfies assumptions («)••• (iv) and if cr and h are bounded Lipschitz con- 
i=i 

tinuous functions, then 

rt rt P rt 

X(i) = X(0)+ / <r(X(s))dW(s) + / b(X(s))ds + V / n i (X(s))dL i (s) (2) 

/ias a unique strong solution with local times Li satisfying I^(-) = / l9D i (X(s)) dLi(s) 

Jo 

Proof of corollary 13.61 

Thanks to proposition I3.4[ V satisfies the assumptions of theorem 13.31 Thus equation 
(J2J) has a unique strong solution X with local time L and reflection direction n. Using 
proposition 13.41 again, there are (non-unique) coefficients Ci(u,s) G [0, 4^] for each normal 
vector in the reflection term to be written as a convex combination : 

n(w,s)= ^ Ci(u,s)ni(K(uj,s)) (3) 

ffl),3X(u,s) 

Let us prove that there exists a measurable choice of the q's : 

The map n (resp. n^(X)) is only defined for (u, s) such that X(cj, s) G dV (resp. X(to>, s) G 
<9T>j). We extend these maps by zero to obtain measurable maps on Q x [0,T] (for an 
arbitrary positive T). Note that equality ([3]) holds for the extended maps too. 

For each (u, s), we define the map f WtS on MP by /u, |S (c) = I S?=i c « n j(X(co>, s))— n(c<j, s)\. For 
a positive integer k, let = {0, |, |, . . . , L^J} P denote the ^-lattice on [0, j^] p endowed 
with lexicographic order. The smaller point in for which f u>s reaches its minimum value 
is 

c (fc) (w, s)=J2 c II (h^Ac^UAc) + I/„,.(c)=/ Ml .(c)I C '>c) 

is a measurable map on x [0,T] and © implies that |/ w ,,(cW(w, s))| < f. Taking 
(coordinate after coordinate) the limsup of the sequence of (c^)*,, we obtain a measurable 
process c^ satisfying (j3J). 

Finally, let Li = lg© i (X(s))cj(s) t£L(s). Lj has bounded variations on each [0, T] since 
the bounded, and it is a local time in the sense of remark 12.11 

Here, strong uniqueness holds for process X and for the reflection term / rij (X(s))dL*(s) 

i=i ^° 

The uniqueness of this term does not imply uniqueness of the Lj's, because the c*'s are not 
unique. ■ 

Proof of proposition 13.41 

Let x G dV. Since Z> = f^ =1 T>i, the set {z s.t. x G <9Dj} is not empty. Since satisfies 
UES{otj) restricted to X>, for each y G Pj (y— x).nj(x) + ^-|y— x| 2 > 0, and consequently : 

Vis.t. &Di3x VyGP (y - x).n i (x) + |y - x| 2 > 

2mm 8 p i3x Q j 
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Summing over i for non-negative q's such that ^g^.^Ci < 4^ we obtain : 
VyGP ( J2 c i n i (x)).(y-x) + - |y - x| 2 > 

which implies that Y^&d ax c « n «( x ) belongs to the set A/x.a of normal vectors on the bound- 
ary of V, for a = /5o mingx'j 9x « j • Thanks to remark this proves the inclusion 

A/^ := < n G S m , n = ^ c i n i( x ) with c i > \ C A4,a 

Let us prove the converse inclusion. We first remark that the smoothness of the boundary 
dT>i at point x implies that for any e > : 

3<^(x) > s.t. {x + z, |z| < 5*(x) and z.n. ; (x) > e|z|} C £>, 
Consequently, for N e = O {z, nj(x).z > e|z|} we have : {x + z, z <E N £ and |z| < 
minj ielcP 

By definition, for each n G A/^f , there exist a n > such that Vy G P (y — x).n + 7^|y — 
x| 2 > 0, hence for z6iV £ and A > small enough : 

A 2 

Az.nH |z| 2 > 

2a n 

For this to hold even with A going to zero, z.n has to be non-negative. So we obtain : 

VneA^ \/e > — n G iV* 

where iV* = {v, Vz G iV e v.z < 0} is the dual cone of the convex cone N e . As proved in 
Fenchel pQ (see also [6]), the dual of a finite intersection of convex cones is the set of all 
limits of linear combinations of their dual cones, in particular : 



For k G N large enough, since — n G A^, there exist unit vectors and non-negative 

k 

numbers q & such that : 



n = ^ Cj,fc n i,fc and for dV l 3 x n.;(x).n ijfc > W 1 - — 
When A; tends to infinity, tends to rij(x), and for fc large enough n^. 1° > f thus : 



1 > n.l° > £ c 4 , fc n 4 , fc .l° > f J] c hk 

9X>3x dX> 4 9x 



Thus the sequences (ci±)k are bounded, which implies the existence of convergent subse- 
quences. Their limits q j00 > satisfy : 

n = Q >00 ni(x) 

d£>i9x 
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This completes the proof of AI~£ C Af^. We already proved that Af^ C AJ~£ a with a = 
Po min,9X) j9 x oij, so we obtain Af£ = Af^ = Af^ a for each x G dV. As a consequence, 
V G UES(P mmi<j< p aj). 

Let us now prove that T> G UNC(/3, S). The Uniform Normal Cone property restricted to 
P holds for Pj, with constant /3j < < 4- That is, for x G Dn <9Pj there exist a unit 
vector 1^ which satisfies l^.rij(y) > \J 1 — f3f for each y G "DPI dVi such that |x — y| < S{. 
If 1^ = n»(x), this implies that |n»(x) — nj(y)| 2 < 2/3?. 

If l x 7^ nj(x), we use the Gram-Schmidt orthogonalization process for a sequence of vectors 
with 1^ and rij(x) as first vectors, then compute n,(x).nj(y) in the resulting orthonormal 
basis (l^,e 2 ,e 3 , . . . , e m ) : 

m(x).m(y) = (^(yM-n^x)) + (n,(y).e 2 ) ^1 - (^(x)) 2 



Note that |rij(y).e 2 | < fa because l^.n^y) > y 1 — (3f and |rij(y)| = 1, thus : 

m(x). ni (y) > (^i-A 2 ) 2 - A 2 = i - 2A 2 

This implies that |n;(x) — rij(y)| 2 < Aj3f. 

So in both cases : |n;(x) — rij(y)| < 2$ as soon as |x — y| < Si for x, y G V fl <9Dj. 

Let us now fix x G <9£> and 5 = mini^p Si/2. We then choose x' G <9I> fl -B(x, 5) such 
that {i s.t. x' G <9£>;} D {z s.t. y G d£>~} for each y G 91? n B(x, S) and we let 1 = 1°,. To 
complete the proof of T> G UNC(P, S), we only have to prove that n.l is uniformly bounded 
from below for n G Afy with y G &D fl -B(x, S). 

We already know that each n G AT^ is a convex sum of elements of the Afy* ■ n = 
^2&d- By c i n i(y)- The coefficients c* are non-negative, their sum is not smaller than 1 
because n is a unit vector, and is not larger than ^ thanks to (iv). So the vector 
n ' = H&D&y c i n i(*') satisfies : 

n.l = Cinj(x').l > fa and |n'-n| < ^ c^n^x')— iii(y)| < 2 ^ c i( 9j < 2 max ^ L 
Consequently : n.l > n'.l — |n' — n| > P — 2 max!<j< p ^ > 0. ■ 

4 Existence of dynamics for globules and linear molecule 
models 

4.1 Globules model 



Let us prove that the globules model satisfies the assumptions in proposition 13.41 The set 
of allowed configurations is : 

a = ( n v a) n ( n n ( n v ^ 

l<i<j<n l<i<n l<i<n 

where = {x G M dn+n , \x,i — Xj\ > %i + Xj] 

V i+ = {x G R dn+n , Xi < r+} Vi_ = {x G M dn+n , £ f > r_} 
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The T>ij have smooth boundaries on A g and the characterization of normal vectors is easy : 
at point x satisfying \x~ xj\ = Xi+Xj > 0, the unique unit inward normal vector n^x) = n 
is given by : 

rii = — — ^ ^— fii = — rij = —4 ^— hi = — (every other component vanishes) 

2(xi + Xj) 2 J 2{xi + Xj) J 2 v 3 1 1 

On the boundary of the half-space T> i+ , at point x such that Xi = r + , the unique unit 
normal vector n has only one non-zero component : hi = — 1. Similarly, is a half- 
space, x belongs to its boundary if Xi = r_, and the unique unit normal vector n at this 
point has hi — 1 as its only non-zero component. These vectors do not depend on x and 
will be denoted by n i+ , rij_ instead of n i+ (x), rij_(x). 



\ I 2 



Proposition 4.1 A a satisfies properties UES = andUNC , . , , 

9 \2r + ny/nJ IV 2 % r\rfi 2 u r\n 6 

Moreover, the vectors normal to the boundary dA g are convex combinations of the vectors 
normal to the boundaries dT>ij, dT> i+ , &D^, that is, for every x in dA g : 

N£ B = < n G S dn+n , n = CijUijix) + c i+ n i+ + c;_n;_ with c ij ,c i+ ,c i _ > 

[ 9X>y3x dT> i+ 3x ai?i_3x 

Proof of proposition 14.11 

We have to check that the assumption of proposition 13.41 are satisfied for the set A g . 

Since the T>i + and are half-spaces, the Uniform Exterior Sphere property holds for 
them with any positive constant (formally «j + = «j_ = +oo). For the same reason, the 
Uniform Normal Cone property holds for T> i+ with any constants /3j + and 5i + , and with 
equal to the normal vector n i+ . This also holds for the sets T>i_ with any and 5j_ , and 
with \ l ~ = n i+ . 

Let us consider x G A g such that x G &Dij, i.e. \xi — Xj\ = Xi + Xj. By definition of rijj(x), 
for y = x — (Hi + Xj)njj(x), one has : 

ry . ry> . ry> . I nn . np . rr* . 

yi = Xi - [Xi + Xj )^-^ = = Xj - [ Xi + xj) 2 ^ i+ ^ = Vj 

y% + Vj = Xi - (xi + Xj)(-\) + Xj - (xi + Xj)(-~) = 2(xi + xj) 

For z G -8(0, Xi + Xj) : 

\{yi+Zi)-(yj+Zj)\-(yi+Zi+yj+Zj) < \zi\ + \zj\-2(xi+Xj) + \zi\ + \zj\ < 2|z|-2(x i +x J -) < 

thus y + z G T>\j . This proves that B(y,Xi + Xj) C V c ip hence N* 3 = AQ* +V The general 
Uniform Exterior Sphere property does not hold, but the property restricted to A g holds 
for T>ij with constant oiij = 2r_ because Xi + Xj > 2r_ for x G A g . 

For x G A g such that x G <9T>ij, let us define 1^ = rijj(x). For another configuration 

y e 92?^ : ' 

F-' n (vi _ Xi — — ^' ~ ^ J i — — — yj ~ | / ^ \2 , / ^ \2 _ ~ x j)-(yj ~ z/j) . i 

x ' ijiyj ~ 2(x i + x J -)"2(y i + J/i) 2(x i + x j )"2(y i + %-) 1 2 J 1 2 J ~ 2(£< + xj)^ + fa) 2 
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Since ( x^ 30 j ^ . {^jji y j ) I <&% Xj I I x<i x j 1 1 \j<i x^ y j ~\~ X j I I X<i X j \ 2 — \/2\xi — x j\\x. — y\ 
and iji + i)j < (x^ + Xj) + \/2\x — y|, this leads to : 

\xi — Xj\ — \/2|x — y| 1 Xj + 



!x ^i(y) > ' : , : , v ' — ^ + r 



2(xi + % + \/2|x - y|) 2 + % + v^lx - y| 

Consequently l^'.ny(y) > y/ 1 — /3 4 2 as soon as |x — y| < (— 7j=r — 1^ • This proves 

that "Djj has the Uniform Normal Cone property restricted to A g with any constant G 

10, 1[ and with <!L- = v^2V- I / 1 « — 1 I > 0. In particular, the property holds for any 

Vv 1 ' / 

with Si, = r„(3yV2. 

For x G dA g , let us construct a unit vector 1° satisfying assumption (iv) in proposition 13.41 
We first have to define clusters of colliding globules : 

^--^x — S.t. j Xj-^ j — X<i I Xj-^ j I Xj-y Xj*2 I — *^3x I «jC 'j'2 5***? I — *^3k I fox SOIX1G 

and their centers of gravity x' = — . . x,-. Let 1° = rr with for each z : 

w xV ; iec x « 

v i =x i -x i and t> j = ( xj 

We need an upper bound on the norm of v. Since x G A g , each Xj is larger than r_ 
and smaller than r + , thus | r ++' r - — x,| < r+ ~ r '" . Moreover, if j belongs to the cluster 
C x (z) around the 2 th globule, the distance between Xj and Xj is at most 2(n — l)r + , thus 
|f j| < 2(ri — l)r + : 

n n 2 

| v |2 = ^ | % |2 + ( r ~ f ^( r+ + r ~ _ x,) 2 < 4n(n - l) 2 r 2 + n— < An 3 r 2 + 

i=l T+ T ~ i=l 

We also need lower bounds on the scalar products of v with normal vectors on the bound- 
aries of the V i+ , V^, T>ij. If Xj = r + , then iij = r+ r _j r r ~2 r+ thus v.n i+ = r -=-. Similarly, if 
Xj = r_, then v.n, = If |x« — Xj| = Xj + Xj then x^ = x^ thus : 

/ \ / /\ Xi—Xj / /\ Xj—Xi 1 r_ /r++r_ v \ 1 r_ /r\f+r_ " \ 

V.rij^XJ — (Xj — ^J- 2 (£. +£j ) T l^Xj — Xj). 2 ^ i +x J ) ~ 2r+~r- \ 2 X V _ 2 r+-r_ V 2 

/ \ Xi—Xj 1 r_ / i o u \ 

— yXi • L 3)-2(xi+$ j ) 2r + -rJ'+ t '- X * X J/ 

— r + (>f _L ^ A _ r - (r++r-) > r__ 
~ 2(r+-r_) ^ X 3> 2(r+-r_) — 2 

because Xj + Xj > 2r_. This proves that assumption (iv) in proposition 13.41 is satisfied with 

Po = — . = =. The inequality /3q > Jl maxfmaxj maxj max, ] /5 i7 ) 

2 v /4n 3 r 2 _ 4r + n-v/n v ,j 

holds as soon as j3y < f3 2 /2. As seen above, the Uniform Normal Cone property restricted 

to A g holds for T>^ with constants (3^ = (3 2 /4 and Sij = j^J=- So, thanks to proposition I3.4[ 
we obtain that A g has the Uniform Exterior Sphere property with constant a^ g = 2r_/9 = 

2r + n^n anc ^ ^ ne Uniform Normal Cone property with constants 5^ = = 2l3 ^ 4 n@ and 



Pa. = V 1 ~ (A/2) 2 = sj 



1 - .7 
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Proof of theorems 12.21 and 12.41 



As seen in the proof of proposition I4.1[ the set A g of allowed globules configurations sat- 
isfies the assumptions of corollary 13.61 If the diffusion coefficients <7j and Oi and the drift 
coefficients 6« and hi are bounded and Lipschitz continuous on A g (for 1 < i < n), then the 
functions 



(71 (Xj 







<ji(x) 








■ 6 x (x) " 






6x(x) 




and b(x) = 









6n(x) 






. 6 n (x) . 



are bounded and Lipschitz continuous as well. Thus, for m = n(d+l) and T> = A g , equation 
([1]) has a unique strong solution X, and there exists a decomposition of the reflection term 
J Q n(s)dL(s) so that it is equal to 




Xi 



X j 



2(X i + X j ) 



{s)dLi 



>).-E 



-dLij{s) - / 
^ Jo 



> + 



KKti 



with 



1 



c, : ,-(s) dLfsl 



L7 



^iW^.Ci+C*) rfL ( S ) and 



1^(8)^+°*- ( s ) ^L(s) for some measurable choice of the Cy, Q+, 



For conve- 



nience, the local time |L y - has been used in equation (£ g ), and denoted by again. 

For an even C 2 function on M^ 1 with bounded derivatives, the function $(x) = J2i<i<j< n fi^i 
Xj) is C 2 with bounded derivatives and the drift function in equation is equal to 
— 7jV$. Thus (5^) has a unique strong solution. Moreover, if Z = J A e~ $(x )c/x < +oo, 
theorem 13.31 implies the time-reversibility of the solution with initial distribution d/x(x) = 
Jl^(x)e-*Wdx. ■ 



4.2 Linear molecule model 

The set of allowed chains configurations is the intersection of 2n + 1 sets : 

n— 1 n— 1 

where A- = {x G M dn+2 , \xi - x i+1 \ > x_} X? l+ = {xG R dn+2 , |x< - x m | < £+} 

P_ = {x6R w L>r_} D + ={xGR w , i + <r + } D = ={x6R w L<i + } 

The boundaries of these sets are smooth, and simple derivation computations give the 
unique unit vector normal to each boundary at each point of dA c (we are not interested 
in other points). Actually, X>_, T> + and T>= are half-spaces, which makes the computations 
and checking of UES and UNC properties very simple. 

• The vector n_ = (0, . . . , 0, 1, 0) is normal to the boundary dT>_ at each point x e <9T>_ 
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Similarly, n + = (0, . . . , 0, 0, —1) is normal to the boundary &D + at each point x £ 
dV+ 

The vector normal to the boundary &D= at each point x £ &D= is n = = (0, . . . , 0, ^= 

At every x £ »4. c fl <9Pi_ (for 1 < i < n— 1), the vector normal to dT>i- is n = nj_(x) 
defined by : 

rii = ^ — -^2±L n . +1 — Xt + 1 x% ^_ _ (other components equal zero) 

x_V3 x_y3 v3 

At every x £ ^4 C D &D{ + , the vector normal to dT>i + is n = n» + (x) defined by : 

Hi = ^ 1 * nj+i = — ^ pi n+ = —= (other components equal zero) 

x + V3 i + V3 v3 



The proof of theorem [23] is similar to the proof of theorems I2.2l and l2.4j. and will be omitted. 
It relies on the following proposition and uses theorem 13.31 and corollary 13.61 to obtain the 
existence, uniqueness, and reversibility of the solution of (S c ). As in the proof of theorem 
12. 2\ the local times are multiplied by a suitable constant to provide a more convenient 
expression. 

So, in order to prove theorem l2.5l we only have to check that A c satisfies the assumptions 
of proposition 13.41 i.e. that the following proposition holds : 



Proposition 4.2 A c satisfies properties UES I I andUNC I Wl 



2rW2V IV 24rV'2 1 M 



Moreover, for each x in dA c , the set J\f^ c of normal vectors is equal to 



nG s dn+2 ,n= ^2 Ci-n»-(x)+ ^2 Ci + n i+ (x) + I gv _(x)c-n- + l dv+ (x)c + n + + I 923= (x)c=n = 

Proof of proposition 14.21 

Let us check that the set A c satisfies the assumptions of proposition 13.41 

The Uniform Exterior Sphere property and the Uniform Normal Cone property hold for 
T>-, T> + and T> = , with any positive a and S and any (3 in [0, 1[, because these sets are 
half-spaces. 

For x £ A c fl dVi_, for a = r ~ 2 , and for y = x — «nj_(x) : 

y . = x .- a — — = x . +1 - a — — y _ = + 

X V *J v V ^ 



so that each z £ £?(0, a) satisfies | (t/i+^j) — (yi+i+^i+i) 
and that |xj — x i+ %\ = x_ implies 



\(yi+Zi)-(y i+ i+z i+1 )\-(y-+z-) < x-(l-^-)+\zi- z i+ x\-x- — — +Z- = -aV3+\zi-z i+ i\+z. 

X — Zi 

We know that \zi — z i+ i \ + z_ < \/3|z| < v^a, thus the above quantity is negative and 
y + z £" Z>j_. Consequently, the Uniform Exterior Sphere property holds for T>i_ with 
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constant = It also holds for V i+ , with any positive constant, because this 

set is convex (the midpoint of two points in V i+ obviously belongs to V i+ ). In order to 
prove that X>j_ has the Uniform Normal Cone property restricted to A c , we fix two chains 
x, y G A c fl 9X>j_ and compute : 

2 Xj - j/j - 1 2|xj - x i+1 \ 2 I 2xi- x i+1 (y { - y i+1 Xj - x i+1 
n i _(x).n i _(y) = 3 . 3 h- = ^ + 



3 x_ y_ 3 3x_ 3 3 x_ \ y_ x_ 

J/j — Xj — Xj+i 



2 

Since \xi — Xj+i| = x_, we obtain : nj_(x).iij_(y) > 1 

3 y- x- 

Then |x_(y; - - y_(x; - x m )| < x_|yj - - X; + x i+ i| + |x_ - y_||xj - x i+ i] leads 
to : 

2 2|y — x| 2|y — x| 

n i _(x).n i _(y) > 1— 7— — - Xi\ + \y i+1 - x i+1 \ + |x_ - > 1 =3 — > 1- 



3t/_ V%_ V3r 



Thus nj_(x).iij_(y) > a/ 1 — Z^ 2 ! as soon as |y — x| < r ^ (1 — a/ 1 — l3f_ ) , which is implied 

by |y — x| < r ~f^ (3f_ . As a consequence, the set £>j_ has the Uniform Normal Cone 
property restricted to A c with any constant G [0, 1[ and the corresponding constant 

Si- = r -=ffi- 

Note that for xje^n 9P i+ : n i+ (x).n i+ (y) = v . — + - 

3 x + y+ 3 

so that the same computation gives that V i+ has the Uniform Normal Cone property 
restricted to A c with any constants (5 i+ G [0, 1[ and 5 i+ = ^-^■f3f + - 

In order to check the compatibility assumption, let us fix a point x G dA c and construct 
a suitable vector v such that 1° = We first define the "middle point" x' of the chain : 

x' = X( n+ i)/2 if the chain contains an odd number of particles, and x' = rr " -/2+ ^" /2+1 if n is 
an even number. We also define x' = x ++ x - , We then construct the v^s in an incremental 
way, starting from the middle and going up and down to both ends : 

• If n is odd, we choose i>( n +i)/2 = 0. 

• In the case of an even n, we choose (v n / 2 , fn/2+1) = (x n /2 — x',x n /2+i — x') if |x n /2 — 
x n /2+i\ < x' and (v n/2 ,v n/2+ i) = (x n/2+ i - x',x n/2 - x') if |x n/2 - x n/2+1 \ > x'. 

In the critical case |x n / 2 — x n / 2 +i| = x', our choice depends on the value of x' : we 
let (v n/ 2,v n / 2+ i) = (x n/2 - x',x n/2+ i - x') if x' < r + and (v B / 2 ,v B / 2+ i) = {x n/2+ i - 
x', x n/2 - x') if x+ = x_ = r+. 

• The other i>j's are chosen incrementally so as to fulfill the same condition : 

- ^i+i = 

Note that \vi\ < (n — l)r + for each i. For our choice of v to be complete, we also define : 

f- — r_ 

• v_ = — , v+ = if x_ < x+ 

2 ' + 2 + 

V 1 

• «_ = _, w+ = - r+ if x_ = x + < r + 

• v_ = — r+, v + = — if x_ = x + = r + 





if 






< x', 


or 


^2+1 1 




if 


\x% 


~ Xj+i 


> x', 


or 


|Xj Xj-)-i | 
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As in the proof of proposition 14. 1[ we need an upper bound on the norm of v : 
|v| 2 = £ H 2 ni^r + f + 1 + ^± < ^ 



8=1 



We also have to prove that the scalar products of v with normal vectors on the boundaries 
are uniformly bounded from below. Note that v.n = V-, v.n + = — v + , v.n = = v+ ^~ , 

v.n i _(x) = — — j=-\Vi - Vi+i) ■= and v.n i+ (x) = — — j=-\v i+ i - Vt) + -=. 

x^yi yS x + yi Vj 

In the case where X- < x+, the choices made on v lead to : 



if x_i 



v.n_ 


r_ 

~~ T 


v.n + 


_ r— 
~ T 








= X + , 



-x i+1 hence v.n. 4 _(x) = > —jf- = ^= 



x + +v + 



> 



We now proceed with the case x_ = x + < r+. Since (a;, — = x^ = x + and t>j — = 
Xi — x i+ i for each i, we have : 

• if X- = r_, v.n = ^ 

3 r_ 

• v.n = = > — != 

, v aL — tL r_ — ^ r_ — x_ + i> + — r + + |r + 

• v.n»_(x) = = — > =r^- = — -= and v.n i+ x = = > — > 



2^3 

In the last case x_ = x + = r + , one has |x, — x i+ i\ = r + and Vi — v i+ i = x i+ i — x^ for each 
i : 

^i + fr+ r+ 



v.rij_(x) = = > — ■= and v.nj+fx) 



\/3 ~ 2y/S y/3 2y/3 

So all these scalar products are larger than As a consequence, assumption (iv) in 

r_ 

proposition 13.41 is satisfied with (3 = -=. Choosing /3j_ = = f3%/4 hence 5j_ = 

r + nyQn 

j^/3q we obtain, thanks to proposition I3.4[ that A c has the Uniform Exterior 



u i+ ~ 64 ^0 

Sphere property with constant a Ac = 7== and the Uniform Normal Cone property with 

Jo 



constants 8 Ac = ^/pr^i and Pa c = \ 1 



24r±n 



2 „:i 
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